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Abstract 

A recently formulated version of the eigenchannel method [R. Szmytkowski, Ann. 
Phys. (N.Y.) 311, 503 (2004)] is apphed to quantum scattering of Schrodinger particles 
from non-local separable potentials. Eigenchannel vectors and negative cotangents of 
eigenphasc-shifts are introduced as eigensolutions to some weighted matrix spectral 
problem, without a necessity of prior construction and diagonalization of the scattering 
matrix. Explicit expressions for the scattering amplitude, the total cross section in 
terms of the eigenchannel vectors and the eigenphase-shifts are derived. An illustrative 
example is provided. 

1 Introduction 

Recently, Szmytkowski pp (see also inspired by works of Garbacz [3] and Har- 
rington and Mautz ^ has presented a general formulation of the eigenchannel method 
for quantum scattering from Hermitian short-range potentials. This formulation, ex- 
ploiting the formalism of integral equations, is alternative to that proposed in 1960's 
by Danos and Greiner 5 , who based on wave equations written in differential forms. 

While various theoretical aspects of the new approach to the eigenchannel method 
have been considered in detail in Ref. Pj, no explicit illustrative applications of the 
formalism have been provided therein. It is the purpose of the present paper to fill 
in this gap by applying the formalism to non-relativistic scattering from non-local 
separable potentials. 

Separable potentials have proved to be useful in many branches of physics such 
as atomic molecular [7] and solid state |5] physics. Still, the most extensive 
use of the separable potentials seems to be made in nuclear physics for describing 
nucleon-nucleon interactions [5]. The utility of separable potentials stems from two 
facts. Firstly, the Lippman-Schwinger |10| equation with a separable potential may 
be solved by employing algebraic methods. Secondly, methods are known which allow 
one to construct separable approximations to arbitrary local potentials [TTj . Here, 
the crucial role plays the method proposed by Ernst, Shakin and Thaler (see also 

lEI)- 

The arrangement of the paper is as follows. Section |21 contains some basic notions 
of quantum theory of potential scattering. In Sectional we focus on scattering from 
separable potentials. In Section 4 we define eigenchannel vectors as eigenvectors of 
some weighted spectral problem and investigate some of their properties. We introduce 



'email; remikOmif. pg.gda.pl 



1 



also eigenphaseshifts using eigenvalues of the same spectral problem. Then, we derive 
formulas for a scattering amplitude in terms of eigenchannel vectors and calculate a 
total cross section and an averaged total cross section. Finally, in Section 5 we consider 
an illustrative example. The paper ends with two appendices. 



2 Non-relativistic quantum scattering from non-local 
potentials 

Consider a Schrodinger particle described by the monochromatic plane wave 

Ur) = (r|fc,) = e"^-'- (2.1) 

with hki denoting its initial momentum, scattered elastically, i.e. |A;,j| = |A;y| = fc [hkf 
is a momentum of the scattered particle) by a non-local Hermitian potential with a 
kernel V{r, r'). 

It was shown by Lippmann and Schwinger (10) that such processes can be described 
by the following integral equation: 

iP{r) = (t>i{r)^ [ d^r' [ d^r" GiE,r,r')Vir' ,r")^{r"), (2.2) 

where G{E, r, r') is the well-known free-particle outgoing Green function and is of the 
form 

ife|r-r'| 

As usual, during considerations concerning scattering processes we tend to find asymp- 
totical behavior of the wave function ip^r). To aim this we firstly need to find asymp- 
totical form of Green function H2.3I) for r — > cxd. After straightforward movements we 
have: 

G{E,ry) .^_e-"=/-, fc, = fc-. (2.4) 



and therefore 



V'(r) ^^e"='-'^ + A(fef ^fe,) — . (2.5) 

r 



The quantity A{kf <— ki) appearing above is called scattering amplitude and is of the 
form 

A{kf ^k,) = --^ f dV / dV'e-''=^ '-V(r',r")^(r"). (2.6) 

In terms of the scattering amplitude the differential cross section for scattering from 
the direction ki into the direction A;^ is 

-^ikf^k,)^\Aikf^k,)\^. (2.7) 

Moreover, one defines the total cross section in the following way 

a{k,)= <f d'kflAikf ^k,)\\ (2.8) 

Jin 

Then averaging the total cross section by all directions of incidence fc,;, one gets the 
so-called averaged total cross section: 

at{E) = d^k i d^kf \A{kf ^ k,)\\ (2.9) 
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3 The special class of non-local separable potentials 

In this section we focus on the special class of Hermitian separable potentials given 
by the following kernel: 

l/(r,r')-5]c.^^;^(r)^;;(r') (V^ : o;^ G M \ {0}) , (3.1) 

where fj, may, in general, denote an arbitrary finite set of countable indices /i^, i.e., 
/i = {/ii, . . . , fXk} and asterisk stands for the complex conjugation. 

Application of Eq. H3.1|l to Eq. H2.2|l leads us to the Lippmann-Schwinger equation 
for the separable potentials: 



V'(r) =</),(r)- Vc.^ / dVG(£;,r,r'K(r') / dV%;(r")V(r" 



). (3.2) 



Similarly, substitution of Eq. H3.1|) to Eq. H2.6|) gives us the scattering amplitude of 
the form 

A{kf^k,) = -^Y.^J dVe-'=/-'t;^(r') / dV't,;(rXr"). (3.3) 

For the sake of convenience, henceforth we shall use the Dirac notation. Then, the 
above formula reads 

A{kf <- fc,) = --^ (3.4) 



The fact that is an element of an arbitrary countable set, implies that we can put 
all scalar products (vfj,\tj}') into a sequence. Therefore, for further simplicity we can 
employ the following notations: 



{v\ip) = 



{Av) I, {'^\v) = {v\^)'' = {{^\vx) {^\v2) . . .) (3.5) 



and 

Q = diag[wi,w2, . . .], (3.6) 

where the dagger denotes the matrix Hermitian conjugation. Note that from Eqs. 
(|3.1() and H3.6|l it is evident that the matrix f2 is invertible. We keep this fact for 
purposes of further considerations. In the light of Eqs. H3.5|l and H3.6|l we may rewrite 
the scattering amplitude (|3.4|l in the following form: 

A{ks^k,)^-^{ks\v)^{v\i,). (3.7) 

In the last step we need to calculate the scalar product To this end, we 

transform Lippmann-Schwinger equation for the separable potentials (|3.2|l to a set of 
linear algebraic equations. Hence, after some elementary movements, we have 

{^vij. + Gi,^LJ^] = (vy\k^), (3.8) 

where 

G.^= / d^r / A^r' vl{r)G{E,ry)v^{r'). (3.9) 
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Finally, application of Eqs. H3.5|l and H3.6f) to Eq. I|3.8f) yields 

(l + GQ)(-i;|V) - (-u|fc,) => {v\^) = {\ + GQ.)-^{v\k,), (3.10) 

where I stands for the identity matrix and G denotes a matrix with the elements 
G,.^. Substitution of Eq. (ITTUIl to Eg. gives the expression for the scattering 

amplitude 

Aikf ^ k,) = -^{kj\v)Q (I + G Q)-' {v\k,). (3.11) 
Since (MN)~^ = N^^M^^ for all invertible matrices M and N, we can conclude that 

A{kf ^ fcO = -^{kf\v) (Q-i + G)"' {v\k,). (3.12) 

4 The eigenchannel method 

Here, we turn to the formulation of the eigenchannel method proposed by Szmytkowski 
0. This author has shown that the eigenchannels provide a powerful mathematical 
tool to the quantum theory of scattering. Moreover, in this approach we do not need 
to construct the scattering matrix to obtain the formula for an averaged total cross 
section. 

In the first step we rewrite the matrix Q.^^ + G as a sum of its Hermitian and 
non-Hermitian parts. Hence, we have 

Q-i + G = A + iB, (4.1) 

where the matrices A and B are defined by the relations: 

A = i [Q-i + G + (Q-i + G)t] = Q-i + i(G + G^) (4.2) 

and 

B = I [n-' + G- (Q-i + G)t] = ^(G - Gt), (4.3) 

respectively. From the definitions (|4.2|1 and H4.3|l it is evident that both matrices A 
and B are Hermitian, i.e., A = A^ and B = B^. Moreover, straightforward calculations 
with the aid of Eqs. (|2.3() and (|3.9() allow us to express their elements as 

[A]., + f d^r f d^r' v:ir rfJ'-- ;'^K ,ir'), (4.4) 

uji. 2TTh^ Jk3 Jk3 r-r' 



sin(fc|r — r' 



[B|.„ = ^/ d-V / i'r'v;{r) ]J _„ " »„('■'). (45) 



Following Jl let us consider the following weighted spectral problem 



AX^{E) = X^{E)BX^{E), (4.6) 

where Xj{E) and XjiE) are, respectively, its eigenvectors and eigenvalues. Through- 
out the rest of the present paper, the eigenvectors X^{E) will be called eigenchannel 
vectors. They constitute a representation of eigenchannels, defined in pp as state vec- 
tors, in subspace spanned by the potential functions v^{r). In fact, the knowledge of 
eigenvectors X^{E) allows us, by some elementary steps, to construct the eigenchan- 
nels. 

By virtue of the fact that both A and B are Hermitian with the aid of positive 
semidefiniteness of the matrix B (for a proof, cf. Appendix 0, we infer that X*{E) = 
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Xj{E) for all 7. Moreover, it is easy to show that eigenvectors X^{E) satisfy the 
weighted orthogonality relation 

XI,{E)BX^{E)=0 (7' 7^ 7)- (4.7) 

Hence, it is natural to assume the following normalization: 

xl,{E)BX^{E)^6y^, (4.8) 

or, using the matrices A and Q^^ + G: 

xl,{E)AX^{E) = X^{E)Sy^, xl,{E){Q-^ + G)X^{E) = [i + X^{E)]Sy^. (4.9) 

Since the eigenvectors X^{E) are solutions of the Hermitian eigenproblem, they should 
obey the following closure relations: 

J2 X^{E)Xl{E)B = I, J2 y^^\E)X^{E)X\{E)k = I, (4.10) 



and 



Henceforth, we shall be assuming that the potential (|3.1|) is such that the above 
relations are satisfied. Therefore, it is possible to express the matrix {Q.^^ + G)~^ in 
terms of the eigenchannel vectors X^{E). Indeed, using Eq. (|4.1HI we have 

(n- + G)-=E^^^.(i^)^^(i^). (4.12) 
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Application of Eq. (|4.12(l to Eq. (|3.12(l yields 



A{kf ^k,)=-^J2 ; + x,{E) {h\'")X,{E)Xl{E){v\k,) (4.13) 
or equivalently 



Because of the very symmetrical form of the scattering amplitude, it is useful to define 
the following functions: 



^-('^^ ^ 'jiSffT.{^\''^)^^^^^) = ^ji^Ak\-)^^^El (4.15) 

hereafter termed eigenchannel harmonics. It follows from their definition that they 
are orthonormal on the unit sphere (cf. Appendix IB|) . i.e., 

d?ky;,{k)y^{k)^5^.^. (4.16) 

After substitution of Eq. (|4.15|) to Eq. 14.1411 , one finds 

A{kf ^ fc,) = -_i^3;,(fc^):y;(fc,). (4.17) 
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Further, it is convenient to express the eigenvalues Xj{E) in terms of so-called eigen- 
phaseshifts 6^ (E) by the relation 

X^{E) = -cot 5-y{E), (4.18) 

which after application to the scattering amplitude H4.17() gives 

A-jr 

A{kf ^ fc,) ^ — Y,^''-^''^ sm8^{E)y^{kf)y;{k,). (4.19) 

7 

As already mentioned the above result was obtained without prior construction of 
the scattering matrix. It is also necessary to emphasize that the method used gives 
formula for the scattering amplitude which has a similar form to that obtained for 
potentials with spherical symmetry V(r) IIQI . For such potentials the functions y^{k) 
reduce to the spherical harmonics Yim{k). Subsequently, combining relation (|4.19|l 
with 1)2. 8|l we have 

"^^i) = -^5]sin2^^(s):^'^(fc,):^^;(fco. (4.20) 

7 

Finally, after substitution of Eq. (|4.19|) to Eq. 12.91) one obtains 

A'jr 

<^t{E) = ^Y.'inH,{E). (4.21) 

7 

Thus, we have arrived at the well-known formula for the averaged total cross section. 



5 Example 

To illustrate the method, let us consider scattering from a pair of identical spheres, 
of radii R. The symmetry of this target allows us to locate the origin of a coordinate 
system in the midpoint of the interval joining the centers of spheres. Thus, we may 
choose the spheres to be centered at points r — ztg. However, due to the assumption 
of non-locality we may simulate this collision process by potential 

V{r,r') = oj [v+{r)vX{r') + v-{r)v*_{r')\ (5.1) 

with 

1 5{\r±e\-R) 

^^^"^=Vi^ ■ ('-'^ 

It should be noticed that potential defined by the Eqs. (|5.1|l and H5.2|l is the special 
case of that proposed recently by de Prunele (see also |H1). As one can notice, for 
further simplicity the strengths of both spheres were taken to be equal and have value 
Lo. Therefore, the matrix Vl defined by Eq. (|3.6|) may be rewritten as 

where I2 is the 2x2 identity matrix. Moreover, straightforward integrations in Eq. 
1)3.9(1 with the aid of the expansion 

I I 1=0 m=~l 

(r> = max{r, r'}, r< =min{r, r'}). (5.4) 
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gives the matrix G in the form 



\MkR)hi+\2kg) hi+\kR) ^ 



(5.5) 



where rj = 2m/fi?, stand for spherical Bessel functions and h^j^\z) stand for 

Hankel functions of the first kind. Then, utihzing Eqs. (|4.2|l and (|4.3|l one finds 

A = 

and 



OJ'^ - vkjoikR)yoikR) ~7^kjl{kR)yo{2kQ) 
-'qkjl{kR)yo{2kQ) - TjkjoikR)yoikR) 

1 Joi2kg)\ 



with yi{z) denoting spherical Neumann functions. Note that, in general 



(5.6) 
(5.7) 



smz 



yoiz) = - 



cosz 



,(+) 



(z) = -i- 



z z 
As one can see the eigenvalue problem (|4.6|) reduces to 

where its eigenvalues X±{E) are 

c^-i - vkjoikR) [yoikR)±jo{kR)yo{2kg)] 



7jkfoikR)[i±Joi2kg)] 



and respective eigenvectors 

X±{E) - {2f^kfo{kR)[l±jo{2kg)]} 
Since for arbitrary k: 

{v+\k) 



{v\k) 



{v-\k) 



-1/2 f 1 
±1 



4T:jo{kR) 



(5.8) 
(5.9) 

(5.10) 
(5.11) 
(5.12) 



and by virtue of Eqs. (|4.15|) and (|5.11|l the eigenchannel harmonics may be expressed 
as 



y±ik) = 



1 



[1±jo(2M]"'^'(c"'"^±g" 



(5.13) 



2\/27r 

and therefore, after substitution of Eqs. H5.10|l and H5.13f) to Eq. (|4.17|) we infer that 
the scattering amplitude is of the form 



A{kf^k,) = -2Tjj^{kR) 



cos{kf ■ q) cos(fci • q) 



u-^ +ij^kjQ{kR)[h^^\kR) + jo{kR)h^t^\2kQ)] 

^ sin(fc/ ■ q) sin(fc, ■ q) 

Lo-^ +if^k]o{kR)[h^^\kR)~ ja{kR)h'^^\2kQ)]_ 

(5.14) 

After application of the above to Eq. H2.8|l and integration over all directions of 
scattering kf we infer 

<J{k^) = 8v^Jo{kR) 

[1 + jo(2fcg)]cos^(fci ■ g) 

{w-i - vkjoikR) [yo{kR)+joikR)yoi2kg)]f + rj^k^jl{kR)[l+jo{2kg)]^ 

[l-jo(2M] sin2(fc,.^) 



{0.-1 - vkjoikR) [yo(fci?) - MkR)yo{2kg)]Y + 7j^k^j^{kR) [1 - M2kg)Y 



(5.15) 



Then averaging Eq. (|5.15|) over all directions of incidence ki and using the fact that 
f d^k, cos^iki ■ q) = 27r[l + jo{2kg)] (5.16) 



'4n 

and 



/ d^fe, sin2(fc, . q) = 2^[1 - joi2kg)], (5.17) 

J4lT 



'4tt 

we obtain the averaged total cross section 

[l+Jo(2M]' 



V-i - r;fcjo(fci?)[yo(fci?) + jo(fci?)yo(2fcg)]}' + r;2fc2j4(fci?)[l + jo(2fc^.)]2 

^ [l-Jo(2fcg)]" 

{w-i - 7ikjo{kR)[yo{kR) - jo{kR)yo{2kg)]}^ + rj^k^j^{kR)[l - jo{2kg)Y 

(5.18) 

Note that the results (|5.17|) and (|5.18|l may be equivalently obtained after application 
of Eqs. lfOn|l and ifO^ . respectively, to Eqs. l|On|l and (|OT|) . 
Acknowledgements I am grateful to R. Szmytkowski for very useful discussions, 
suggestions and commenting on the manuscript. Discussions with P. Horodecki are 
also acknowledged. 

A Positive semidefiniteness of the matrix B 

The proof is due to R. Szmytkowski 17 . Below we show that the matrix B given by 
Eq. (|4.3|) is such that, the inequality 

Xl{E)BX^{E)>0 (A.l) 

is satisfied. Since 

/ d^fce'^-(— ')=47r^^^ifcfl (A.2) 
Ji^ \r - r'\ 

we may rewrite Eq. (|4.5() in the form 

[Bl^'M = /Si / d^fe / d=^r<(r)e'^- / dVx;,(r')e-"=-', (A.3) 
which after application to Eq. (|A.ip yields 

XliE)BX,iE) = ^£d2fc ^X;,iE)l^^d'rvtir)e^ 
Obviously, the above statement finishes the proof. 



k-r 



2 



> 0. (A.4) 



B Proof of orthonormality relation ( 14.161) 

Substituting Eq. (|n5|l to Eq. l|^?TH|l and reorganizing terms one finds 

(f d?ky;,{k)y^{k) 

J Air 



X 



mk 
8^ 



Y,X;,,{E) [ d?rv:{r) [ d^r' v^{r') <f d'ke"'<^-^">X^^{E). 

JB3 JR3 JiTT 

(B.l) 
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In virtue of Eq. lfO|) . Eq. l|n?T|l may be rewritten in the form 




(B.2) 



As one can notice with the aid of Eq. (|4.5|l . the expression in square brackets is 
proportional to an element of the matrix B. Hence, we arrive at 



which after comparison with Eq. (|4.8|) completes the proof. 
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